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Abstract 

The soft modes in the two- and eight-direction order-disorder ferroelectrics 
are calculated under the mean-field approximation. We find that the conven- 
tional method of the pseudospin model to calculate the soft-mode frequency 
is incorrect, and present a valid modified method. It is demonstrated that the 
conventional method does not show the soft-mode frequency going to zero 
at a critical temperature in the presence of random internal fields, while the 
frequency calculated by our modified method goes to zero in random fields 
at a critical temperature. In the eight-direction ferroelectrics, the soft-mode 
frequency decreases to zero at the first-order phase transition temperature 
though the symmetry has been destroyed at high temperatures under an ex- 
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ternal field. The promotion effect of random fields on the phase transition is 

testified by the calculation results of the soft modes in the paraelectric phase. 
PACS: 77.80.-e 64.60.-i 63.70.+h 77.80.Bh 
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I. INTRODUCTION 



A pseudospin model is usually used to interpret the phase-transition phenomena of order- 
disorder ferroelectrics such as KH 2 P04 in terms of the microscopic atomic interactions.0 The 
model describes the ferroelectric phenomena as the motion of active ions in the double-well 
type potential field. The dynamic properties are generally considered by calculating the 
frequency of the ferroelectric soft mode. It yields satisfactory results for the second-order 
ferroelectric phase transition: when the temperature increases, the frequency of the soft 
mode decreases below the critical temperature T c and increases above T c . At the critical 
temperature T c , the frequency is equal to zero, i.e., the ion vibration is frozen. In the 
previous works,Bi the conventional equation to determine the frequency of the soft mode is 
expressed as 



d(Si) t 



-i([Si,H]) t , (1) 



dt 

where S« is the pseudospin of the i-th ion, H is the Hamiltonian of the system, and t is 
the time. (. . .)t denotes the temporal thermal averaging. The equation (1) is obtained by 
making thermal averaging on the quantum mechanics equation 

f = -M (2) 

To obtain Eq. (1) from Eq. (2), an assumption implied here is 

d(Si) t = (dSi) t , (3) 

i.e., 

d{Si) t = ^2p a dS ia , (4) 

a 

where S ia is the spin of the a-th eigen state, and p a is the statistic weight. However, starting 
from the original definition of (Sj)t, i.e., 

(Si)t = Pa^ia, (5) 

a 



the correct equation should be 

d(Si) t = Y^[padS ia + S ia dp a ). (6) 

a 

Therefore, one purpose of this article is to investigate the soft mode of the order- disorder 
ferroelectrics by using Eq. (6) and compare with the results of the conventional method. 

Rather different from the normal ferroelectrics, relaxor ferroelectrics (relaxors) experi- 
ence no macroscopic phase transition at zero external field.i'i However, a strong external 
electric field can induce a first-order ferroelectric phase transition in the relaxors such as 
Pb(Mgi/3Nb2/3)03 (PMN).H The phase transition and other special characteristics (such 
as the strong frequency dispersion and the history-dependent polarization behaviors) of re- 
laxors are widely believed to be caused by the random interactions and the random electric 
fields in the system. Based on the experimental observationS that the polar active ions 
in PMN shift along eight {lll}-equivalent directions, an eight-direction order-disorder fer- 
roelectric modeli was presented to explain the induced first-order phase transition. The 
dynamic properties, however, are not calculated in that workjll So, another purpose of this 
article is to examine the dynamic properties of the phase transition in the eight-direction 
order-disorder ferroelectric system. 



II. METHODS 



In the eigen state space of the operator S z , the Hamiltonian of the two-direction pseu- 
dospin model can be equivalently expressed as 



H 



-E Po § 



n 
2 



(7) 



where E is the electric field strength, and po is the magnitude of the dipole moment when 
an ion locates in a certain potential well. Q is the tunneling frequency between the two 
potential wells, whose magnitude is equal to the energy difference between the asymmetric 
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and the symmetric zero field. The dipole moment p is associated with S z 

through the relation 

P = 2p S z . (8) 

Under the mean field approximation, the interactions upon a certain ion may be represented 
by an equivalent field, i.e., 

(p) 

Ep = J h E ext p + E rand p , (9) 

Po 

where (p) is the thermal average value of the dipole moment, and J is the ferroelectric 
coupling energy between polar ions. E ext is the applied external electric field, and -E ran d is 
the internal random fields in the system. -E ran d is necessary to explaining the phase transi- 
tions in incipient ferroelectrics such as KTa03:Li,Nb,Na; SrTi03:Ca or PbTe:Ge with dipole 
impurities (active ions with random sites and orientations) .0 It comes from the direct inter- 
actions of electric dipoles, fields created by point charge defects and aforementioned defects, 
etc. Random fields are also important to interpreting the special properties of ralaxors such 
as PMN.EH§ The pinning effect in PMN was believed to be induced by the random fields 
of impurities in the system.i In the process of the thermally activated flips of the local 
spontaneous polarization, the random interactions between polarization microregions are 
essential in producing the ralaxor characteristics such as the diffused phase transition and 
the frequency dispersion.! It is also revealed that the charged chemical defects and nano- 
domain textures originating from atomic ordering also give important contributions to the 
random fields.lli The real occurrence of long-range ordered phase depends on the competi- 
tion of the constant-sign ferroelectric coupling and the alternating-sign random fields. The 
amplitude of random fields is mainly determined by the impurity and the microstructure 
imhomogeneity in the system. If there are more charge impurities or the microstructure is 
more inhomogeneous, the random fields are stronger. Random fields coming from different 
sources have different distribution,!'^ but have similar influences. In this article, £7 ran d is 
assumed to have a Gaussian distribution with a width a P as 



P(-Erand) 



2-ku, 



exp 



-E 2 

rand 



(10) 



The static properties of the system can be easily determined from Eqs. (7-10), which have 
been discussed extensively previously. 00 

To investigate the dynamic properties of the system, the fluctuation of the electric field 



dE(t) = dE ■ exp(iu;t) 



(11) 



is considered. The eigen energy and the eigen polarization (dipole moment) of the Hamil- 
tonian in Eq. (7) are respectively expressed as 



E ± = ±J(E Po ) 2 + (n/2y 



(12) 



and 



P± 



Epo_ 



-Po- 



(13) 



The variation of p± induced by dE(t) can be calculated in a linear perturbation method of 
the quantum mechanics, which yields 

n 2 E±p dE(t) 



dp±(t) 



Po- 



{Epo) 2 AE 2 ± -{huj) 2 
According to the conventional method [Eq. (4)], d(p) t reads as 



(14) 



d(p)t 



dp + (t) tanh 



■ p(E ra , n d)dE Tan( i, 



(15) 



while in the modified method, d(p)t reads as 



d(p)t 



dp + (t) tanh 



E + 
kv.T 



+ P+ 



1 - tanh 2 



E + 
kv.T 



Ep 2 dE(t) 
k B TE + 



P(-Erand)d-Erand- (16) 



Thus the frequency of the soft mode, u, can be determined from the self-consistent equation 

, d (p)t 



p dE(t) = J- 



Po 



(17) 



The dynamic properties of the pseudospin model have been investigated in details by using 
the conventional method when there is no random field.0! 



For the eight-direction order-disorder ferroelectric model, there are eight potential wells 
along {lll}-equivalent directions, and the Hamiltonian is assumed to be the direct sum of 
four groups of double- well matrixS, i.e., 



H 















-\Ep § 











n 

2 








lEp 















■\Epo f 

f \Ep, 

-±Ep 





n 

2 











n 

2 

jEpo 



(18) 



where the factor 1/3 is the cosine of the angle between the dipole moment and the external 
field when they are not parallel to each other. Other equations similar to the two-direction 
cases can also be obtained, which are not listed here in order to keep the article concise. 



III. RESULTS AND DISCUSSIONS 

A. Two-direction pseudospin model 

In order to demonstrate the different effects of the conventional and the modified meth- 
ods, the soft-mode frequency of the two-direction pseudospin model corresponding to Q = J 
is shown in Fig. 1 when there are no external field and random fields. The dashed line 
represents the calculation results by using the conventional method, and the solid line repre- 
sents those calculated by using the modified method. We can see that both the conventional 
method and the modified method yield the same frequencies of the paraelectric phases 
(T > T c ). The reason is that the polarization and z-spin of the eigen states in the para- 
electric phase are zero [p± = while E = 0. See Eq. (13)], so the Eq. (4) is equivalent to 
Eq. (6) in this case. For the ferroelectric phase, {Tiuj) 2 in the conventional method is pro- 
portional to the polarization of the system;0 (hu) 2 in the modified method is smaller than 
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that in the conventional method. The feature can be understood as follows: under the same 
drive dE(t), the "displacement" d(Si) t calculated by Eq. (4) is smaller than that calculated 
by Eq. (6), i.e., the "elastic" coefficient of the conventional method is larger than that of 
the modified method, so the vibration frequency in the conventional method is higher than 
that in the modified method. It should be noted that the frequency decreases to zero at the 
critical temperature, which is just a unique characteristic of the soft mode theory. 

Figure 2 shows the polarization and the soft-mode frequency of the system under a weak 
external field E ext = 0.005 J/po when Q = J. It can be seen that the polarization increases 
continuously when the temperature decreases. The soft-mode frequency does not decrease 
to the zero value but reach a nonzero minimum. The minimal value of the frequency in 
the modified method is closer to zero than that in the conventional method. The frequency 
calculated in the modified method reaches its minimum at = 0.95 J/k^, which is a little 
larger than the temperature at which the polarization changes most quickly, T m = 0.91J/fcB- 

The random fields existing in the system have great influence on the properties of disor- 
dered ferroelectrics.0 They can make the critical temperature decrease, and even inhibit the 
phase transition. The polarization is calculated for the two-direction model with distribu- 
tion width of random fields a e = and 0.3 J/po, respectively, and the results are presented 
in Fig. 3(a). It is demonstrated that the critical temperature T c decreases from 0.911 to 
0.802 J/k-Q when the random field width increases from to 0.3 J/po- The figure 3(b) 
shows the calculated results of the soft-mode frequency when a e = 0.3J/p and Q = J. We 
observe an extraordinary difference between the results calculated by the two methods: the 
frequency in the conventional method keeps nonzero while the frequency in the modified 
method decreases to zero at the critical temperature. From the curve of the polarization in 
Fig. 3(a) we know that the phase transition is of the second-order type and the symmetry of 
the system changes at the critical temperature. A zero frequency at the critical temperature 
is required by the soft mode theory. Thus, Eqs. (1) and (4) used in the past to calculate 
the frequency are incorrect, while Eq. (6) produces the physically reasonable results. 
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B. Eight-direction order-disorder ferroelectric model 

Under the two-body couplings, there is only the second-order phase transition in the 
two-direction ferroelectric system, while both the second-order and the first-order phase 
transitions can exist in the eight-direction ferroelectrics.0'0 The polarization and the soft- 
mode frequency in both cooling and heating processes are shown in Fig. 4 for the eight- 
direction model when Q = 0.63 J. Fig. 4(a) clearly demonstrates the jump of the polarization 
and the difference of the phase transition temperature between the cooling and heating 
processes, which indicates that the phase transition is of the first-order type. In a cooling 
process, the soft mode frequency decreases to zero at the phase transition temperature (over- 
cooled temperature) and suddenly jumps to the nonzero frequency value of the ferroelectric 
phase. In a heating process, the cases are similar except that the phase transition occurs at 
a higher temperature (overheated temperature). 

For the relaxor ferroelectrics such as PMN, an external field with proper strength will 
induce a first-order phase transition.i'@ The second-order phase transition occurs in stronger 
field ranges, and the phase transition vanishes when the field increases further. The curves of 
the soft-mode frequency at different external fields are depicted in Fig. 5 when Q = J. At the 
external fields E = 0.075 and 0.095 J/po, the phase transition is of the first-order type (see 
Ref. 12). There is a jump of frequency at the phase-transition temperature, and a down-pulse 
reaching the zero value exists above the phase transition temperature. When an external field 
is applied, the symmetry of the system has already been lowered at high temperatures, and 
the high and low temperature phases can not be distinguished from the symmetry changing. 
However, the zero soft-mode frequency at the critical temperature clearly indicates the 
existence of the phase transition. For E = 0.115J/po ; the phase transition is of the second- 
order type, and a down-peak of the frequency is easy to be identified. It is not a "real" phase 
transition in the sense that the soft-mode frequency does not decrease to zero. The phase 
transition can only be approximately defined from the rapid change of the polarization or 
from the minimum of the soft-mode frequency. For a strong field E = 0.l5J/po, the soft- 
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mode frequency varies rather smoothly, and the softness of frequency is indefinite. In this 
case, it is difficult to identify any phase transition; in other words, the transition does not 
exist. 

Rather different from the common sense, the random fields with proper distribution width 
will promote the spontaneous phase transition in the eight-direction ferroelectric system.0 
In Fig. 6, the soft-mode frequencies of the paraelectric phase at different temperatures when 
Q = J are given as functions of the width of random fields. With increasing width of the 
random fields, the frequency decreases first, and then increases. The minimal frequency is 
reached at a nonzero width. The promotion effect of random fields on the phase transition is 
demonstrated most clearly in the case of T = 0.085 J/k^: the random fields with distribution 
width 0.276J/po < cr c < 0.436J/po induce the spontaneous appearance of the ferroelectric 
phase. 

The original motivation of the eight-direction order- disorder ferroelectric model is to 
predict the phase transition in relaxors, but not to explore the dielectric mechanism. It 
is known that relaxors have a rather broad and smooth spectrum of the relaxation times, 
and the relaxation time shows a progressive slowing down upon cooling, which results in a 
transition into a glassy state.i'ill A distribution of relaxation time instead of a single fre- 
quency of the soft mode should be introduced for fully describing the dielectric properties of 
relaxors such as the frequency dispersion. The Monte Carlo simulation on the two-direction 
pseudospin model has been done to explain the dielectric behaviors of relaxors.10 Thus 
further works of simulation on the eight-direction model may be helpful for understanding 
the characteristics of relaxors. 

IV. SUMMARIES 

In summary, the soft modes in the two-direction and eight-direction order- disorder ferro- 
electric systems are syudied in this article. We indicate that the conventional method yields 
incorrect soft-mode frequency of the pseudospin model, and present the valid modified equa- 
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tion. The soft-mode frequency of the ferroelectric phase calculated by the modified method 
is smaller than the conventional values, while the critical temperature does not vary. When 
there exists random electric fields, the frequency calculated by the conventional method does 
not decrease to zero at the critical temperature, which is physically unreasonable. In con- 
trast, the modified method yields the reasonable results. In the eight-direction ferroelectrics, 
the soft mode of paraelectric phase decreases to zero at the over-cooled temperature. The 
calculation of the soft modes under external electric fields shows that the frequency de- 
creases to zero at the first-order phase transition temperature although the symmetry of 
the system has been destroyed at high temperatures. We find that the soft-mode frequency 
decreases first and then increases with increasing distribution width of the random fields, 
which testifies the promotion effect of random fields on the spontaneous appearance of the 
ferroelectric phase. 
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FIGURES 

FIG. 1. Soft-mode frequency of the two-direction pseudospin model as function of temperature 
when f2 = J. The dashed and the solid lines represent the results in the conventional and the 
modified methods, respectively. The temperatures is measured in unit of J/ks, and hu in unit of 
J. 

FIG. 2. Temperature dependence of the polarization and the soft-mode frequency of the 
two-direction ferroelectrics at an external field E ext = 0.005J/po when 17 = J. The polariza- 
tion (p) is measured in unit of po. The dashed and the solid lines represent the results in the 
conventional and the modified methods, respectively. 

FIG. 3. (a) The polarization of the two-direction ferroelectrics when a e = (dot-dashed) and 
0.3 J/po (solid), respectively, (b) Soft-mode frequency of the two-direction ferroelectrics as function 
of temperature when Q = J and a c = 0.3 J /po, where the dashed and the solid lines represent the 
results in the conventional and the modified methods, respectively. 

FIG. 4. Temperature dependence of (a) the polarization and (b) the soft-mode frequency of 
the eight-direction ferroelectris when Q = 0.63J. The solid and dashed lines represent the results 
in a cooling process and a heating process, respectively. 

FIG. 5. Soft-mode frequency of the eight-direction ferroelectrics as function of temperature at 
a cooling process when Q, = J. Curves 1-4 correspond to E cxt =0.075, 0.095, 0.115, and 0.15 J, 
respectively. 

FIG. 6. Soft-mode frequency of the eight-direction ferroelectrics as function of random field 
distribution width when VL = J . The solid, dashed, and dot-dashed lines correspond T=0.1, 0.09, 
and 0.08 J/k^, respectively. 
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